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Calculating the nuclear mass at finite angular momenta 
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Mean field methods to calculate the nuclear mass are extended into the high-spin regime to cal¬ 
culate the nuclear binding energy as a function of proton number, neutron number and angular 
momentum. Comparing the trend as a function of mass number for a selection of high-spin states, 
a similar agreement between theory and experiment is obtained as for ground state masses. 
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A fundamental property of nuclei is their mass or 
equivalently, their binding energies, B. The variation of 
the nuclear mass with proton and neutron number will 
reveal the shell effects which are closely related to the 
magic numbers and the extra binding associated with 
these numbers. It will also give some general idea about 
which regions of nuclei are deformed and if some specific 
particle numbers gives rise to extra binding for deformed 
nuclear shapes. In recent years, it has become possible 
to study a large number of nuclei up to very high angu¬ 
lar momenta. A natural extension is then to study the 
variation of the total nuclear energy as a function of the 
angular momentum, I, i.e. to extend the investigations 
of the binding energy to three dimensions, B{Z, N, I). 

A first attempt to study the experimental shell effects 
at high spin was carried out in Ref. Q where the energies 
of high-spin states in Z = 50—82 nuclei were plotted rela¬ 
tive to a somewhat schematic rotating liquid drop energy. 
In the present study we will instead investigate how well 
the high-spin states in a few selected nuclei are repro¬ 
duced by state-of-the-art macroscopic-microscopic calcu¬ 
lations. A big simplification at high spin is that pairing 
correlations are negligible which should make our calcu¬ 
lations more reliable. We will thus consider a limited 
number of nuclei whose level schemes are known up into 
the unpaired regime and which have been successfully 
interpreted in calculations. 

In the present study, we will start from the finite range 
version of the liquid drop model because for this model 
systematic fits to the masses of all nuclei with N,Z > 8 
have been performed Q. It turns out however that this 
model has some problems at low mass numbers where the 
energy becomes unstable with respect to high multipole 
deformations 0. Therefore, we will also consider a recent 
version of a classical liquid drop formula, the LSD model 
0 where an curvature term has been added to the 
classical Myers-Swiatecki expressions 0. 

In the macroscopic-microscopic approach, the total nu¬ 
clear energy is obtained as 

Ei.nXdef) + + i?shen(de/, I) 
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in the high-spin limit when pairing is ignored. In the 
formula, it has been specified which terms depend on de¬ 
formation {def) and angular momentum, I, respectively. 
The two first terms correspond to the rotating liquid drop 


energy and the third term is the shell energy. 

In order to calculate the shell effects at high spin, we 
will rely on cranked Nilsson-Strutinsky (CNS) formalism 
with the modified oscillator potential because it is only 
in this model HO that systematic high-spin calculations 
have been carried out for nuclei in essentially all mass 
regions. In order to keep the number of parameters as 
small as possible, we have used the so-called A = 110 
parameters for all nuclei. These parameters have been 
optimized for nuclei with A — 100 — 150 but should be 
approximately applicable for all mass numbers. 

As one alternative for the static liquid drop energy. 
Eld. in Eq. iP), we consider the finite range liquid drop 
model (FRLDM) Q which is more consistent 0 than the 
finite range droplet model (FRDM) in the description of 
fission and thus probably more reliable for strongly de¬ 
formed shapes. The FRLDM parameters have been fitted 
to reproduce ground state masses for 1654 nuclei and 28 
fission-barrier heights. The resulting root mean square 
error after experimental errors have been compensated 
for is (Tth = 0.779 MeV 0. In the present study, this 
model can not be applied directly because pairing is im¬ 
portant at the ground state while we are interested in 
unpaired high-spin states. Indeed, in the fit, an aver¬ 
age pairing energy is included in the macroscopic liquid 
drop energy and must be removed. Furthermore, a zero- 
point energy for vibrations in the elongation direction is 
also included which appears inconsistent for high-spins 
where the special symmetry favoring axial shapes is bro¬ 
ken. One might consider to include the full quadrupole 
zero-point energy term but this would be very difficult 
in practice. Thus, a new fit to the same masses has been 
performed where the full pairing interaction has been in¬ 
cluded in the microscopic terms and where no zero-point 
energy is included. This new fit gives a very similar mean 
square error, namely ati = 0.778 MeV, corresponding to 
a ‘standard’ root mean square error of 0.783 MeV. The 
biggest difference between the two fits is the constant 
term which increases from 2.6 to 7.2 MeV but there are 
also minor differences in all other terms. 

The second alternative for the static liquid drop en¬ 
ergy, the LSD model 0, has been fitted to the same 
shell corrections as the FRLDM for Z > 29 and N > 29 
but for lighter nuclei semi-empirical shell corrections are 
used instead. Therefore, as illustrated for nuclei along 
/3-stability in Fig. QJ the two models give very similar 
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Figure 1: Difference in energy between the FRLDM and the 
LSD model for nuclei along the beta-stability line at spin zero 
and spherical shape. 


energies for mass numbers ^ « 50 — 250 but they be¬ 
come rather different for smaller mass numbers. The 
LSD model has been fitted to 2766 masses with a root 
mean square error of 0.698 MeV which decreases to 0.610 
MeV when only 1654 masses are considered. In order 
to remove the zero-point energy and the average pairing 
when using the LSD model without refitting the param¬ 
eters we make an estimate based on the two fits for the 
FRLDM. We calculate the energy difference for a spheri¬ 
cal nuclei between the two fits and add this contribution 
to the LSD model. This way we introduce a constant 
shift for each nucleus but we do not change the original 
deformation dependence in the model. 

With the static liquid-drop energy fixed, it only re¬ 
mains to fix the parameters of the rigid moment of iner¬ 
tia which enters in the second term of Eq. O- It can be 
calculated as @ 

Jrig. = ‘^MrlA‘^^^5j{£2,j,e4) + 4:Ma‘^ (2) 

where the radius for a spherical nucleus is parameter¬ 
ized as and where the diffuseness is introduced 

in the form of a Yukawa folding function with range 
a. A density distribution which changes from 10% to 
90% of its central value in a distance of 2.4 fm is de¬ 
scribed by a « 0.75 fm. M is the nuclear mass which 
varies as A. The first expression on the right in Eq. 
0 is the rigid body moment of inertia for a spherical 
nucleus with a sharp boundary. The deformation correc¬ 
tion, Sj-(£ 2 , 7 ,£ 4 ), which equals one for spherical shape, 
is larger than one for shapes mainly relevant for rapidly 
rotating nuclei, i.e. for rotation around the smaller axis 
(0° < 7 < 60°). Note that the diffuseness correction (the 
last term in Eq. is independent of deformation. 

The moment of inertia parameters, rg and a are ob¬ 
tained by a fit to experimentally determined nuclear 
charge density distributions. We fit the root-mean-square 
value of the radius (r^}^^^(£ 2 ,ro,a), calculated with a 
Yukawa folding function, to the values given in Ref. 0 
for 116 nuclei with A > 16. The quadrupole deforma¬ 
tions of the ground-states are taken from Ref. Q. The 
result of the fit is rg — 1.1599 fm and a = 0.5984 fm 
and the standard deviation of the errors are s = 0.0454 
fm. Since the root-mean-square radius as well as the mo¬ 
ment of inertia involves integrals over (x^), {y^) and {z^), 
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Figure 2: Experimental (top panel) and theoretical (middle 
panel) microscopic energies and their difference (lower panel) 
for some well-established high-spin bands in A = 109 — 156 
with the level scheme of ^°®Sb extended down to 7 = 0. The 
LSD model and a diffuse surface (ro = 1.16 fm, a=0.6 fm) has 
been used when calculating the first and second terms, respec¬ 
tively, of Eq. Q. Calculated states where all valence nucleons 
have their spin vectors aligned with the axis of rotation (ter¬ 
minating states) are encircled. 


this should give accurate values for the moment of inertia 
even though there might be a larger uncertainty in the 
two values, rg and a. 

Fig. [21 illustrates the total energy for a few nuclei as a 
function of spin. It is drawn using the LSD model for the 
static liquid drop energies and radius constants ro=1.16 
fm and a=0.6 fm obtained from the fit described above. 
For each nucleus the energy of the corresponding rotating 
liquid drop has been subtracted. This is thus a straight¬ 
forward generalization of ‘standard’ mass plots, see e.g. 
Figs. 1 and 2 of Ref. 0, which gives us a convenient scal¬ 
ing when comparing the microscopic energy for different 
nuclei or when comparing theory and experiment. 

The yrast line for ^°®Sb in Fig. O is drawn starting 
from low spins where the discrepancies between between 
calculations and experiment should be an approximate 
measure of the pairing energy which is not included in 
the calculations. The high-spin bands in are 

both observed to termination at 7 = 41.5 and I — 45, 
respectively. They are some of the first and best devel- 
OTed smooth terminating bands in the A = 110 region 
Q. Relative to the ^°°Sn core, they have two proton 
holes in the gg /2 shell and one proton and two and three 
neutrons, respectively, in the high-j hii /2 orbitals. The 
band in ^^^I, referred to as band 1 in Ref. , has the 
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FRLDM, 1 ;,= 1.16 fm, a = 0.60 fm 
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Figure 3: Difference in total energy between theory and ex¬ 
periment for 10 nuclei with masses in the range A — 20 — 194. 
In the lower panel, the same parameters as in Fig.|2|have been 
used while in the upper panel, the LSD model for the static 
liquid drop energy has been replaced by the FRLDM and, in 
the middle panel the diffuse surface when calculating the rigid 
moment of inertia has been replaced by a sharp surface. 


same number of holes and thus more valence particles 
which means that it has a higher terminating spin value 
at = 54.S’*" and it is observed up to / = 50.5. The 
i56Dy band 0 shown in Fig. [21 is yrast at the highest 
spin values. It is interpreted as having four proton holes 
in the Z = 64 core and is tentatively observed to its ter¬ 
mination at / = 62. The band drawn for ^®^Dy is the first 
superdeformed band identified 0 in the A = 150 region 
and the only such band whose absolute excitation energy 
is known [3| • Its interpretation in CNS calculations is 
discussed in Ref. 0 - 

In order to extend the comparison of Fig.j^lto a larger 
mass region, we compare calculated and experimental en¬ 
ergies for additional nuclei in Fig. 0 where different ver¬ 
sions of the rotating liquid drop formula are used. Before 
discussing the details of this figure, we will briefiy de¬ 


scribe the rotational bands which in addition to the ones 
previously discussed have been included in the compar¬ 
ison. The band observed to terminate at 8 + in foNeio 
[l^ is built with two valence protons and two valence 
neutrons in the shell. Cranked Nilsson-Strutinsky 
calculations for this band have been discussed e.g. in Ref. 
0 and reviewed in 0|. The band in ^®Ar is the so called 
superdeformed band which is interpreted in CNS as well 
as shell model calculations as having identical proton and 
neutron configurations with two particles in i'j /2 and two 
holes in P 3/2 0 - The highest spin band shown for ®®Cu 
0 has been referred to as superdeformed with four holes 
in the ®®Ni core and consequently seven particles in open 
shells with three of them in high-j g 9/2 orbitals. The 
other band (with both signatures shown) has three holes 
and two of the six particles in g 9/2 orbitals. The two 
bands in ®^Zn have one hole in the core and two or three 
of the seven valence particles in the g 9/2 shell [ 2 ^. For 
^®®Dy, the two negative parity bands observed to highest 
spin, = 52“, 53“ are included in addition to the band 
shown in Fig.O Of the three superdeformed bands linked 
to the normal-deformed states in the A = 190 region, the 
band in 0 has been included in our compari¬ 

son because it is observed to higher spins than the other 
bands. According to the present calculations, it has the 
high-j 7 r(ii 3 / 2 )'^i 2 (ji 5 / 2 )"^ configuration in agreement with 
standard interpretations, see e.g. j^. 

With the same parameters used in Fig. O and in the 
lower panel of Fig. (31 those data points that are repeated 
are identical. These parameters describe the data with 
good accuracy. Typical errors are in the range of ±1 
MeV. Note especially that from the limited data set con¬ 
sidered here, the general trend as a function of mass num¬ 
ber appear correct. This is contrary to the results in the 
upper panel where in the static liquid drop energy, the 
LSD model has been replaced by the FRLDM. This re¬ 
sults in systematic differences for light nuclei where the 
calculated energies are lower than the experimental en¬ 
ergies. Thus, the FRLDM predicts similar results as the 
LSD model for heavy nuclei but substantially lower total 
energies for ^°Ne, ^®Ar and ®®Cu. One reason for this is 
that the two models are fitted to different shell correc¬ 
tions for Z <29 and N <29 (as seen in Fig. QJ. Another 
reason is that for light nuclei the two models have differ¬ 
ent deformation dependencies, with the FRLDM being 
softer in the £2 and £4 directions. Therefore for light 
nuclei the FRLDM generally predicts somewhat larger 
£2 deformations and in some cases so large £4 defor¬ 
mations that the corresponding shapes appear unreal¬ 
istic. Consider for example the terminating state in ^®Ar 
which is predicted to have a deformation of £2 = 0.415 
and £4 = 0.188 using the FRLDM and £2 = 0.356 and 
£4 = 0.065 in the LSD model. Similar results are ob¬ 
tained for the terminating state in ^°Ne which is pre¬ 
dicted to have a deformation of £2 = 0.120 and £4 = 0.210 
using the FRLDM and £2 = 0.092 and £4 = 0.033 in the 
LSD model. 

The only difference between the two lower panels in 
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Fig.El is that a sharp surface is used when calculating 
the rigid body moment of inertia in the middle panel 
while a diffuse surface is used in the lower panel. The 
radius in the sharp surface case is chosen as the value 
used to calculate the Coulomb energy in the LSD for¬ 
mula, To = 1.21725 fm corresponding to a rigid moment 
of inertia of a sphere which is about the same same as 
in the diffuse surface case for A w 140. In general, the 
differences between the two lower panels are not so big 
for heavy nuclei but become more pronounced for light 
nuclei. For example, the ‘diffuse moment of inertia’ for 
spherical shape is around 13 % bigger for A = 36 and 
more than 20 % bigger for A = 20. For ^°Ne and ^®Ar, 
this results in an energy difference between 2.5 and 3 MeV 
for the terminating 1 = 8 and and 7 = 16 states, respec¬ 
tively. For large deformations however, these differences 
become smaller because the diffuseness correction to the 
moment of inertia is independent of deformation (Eq. |2l . 
Thus, for the strongly deformed 16+ state in ^®Ar, the 
difference is only about 1.5 MeV in the full calculation 
while for the less deformed 8+ state of ^°Ne its around 
3 MeV. Differences up to about 1 MeV are then seen in 
the bands for the A = 60 nuclei. An interesting effect 
is seen in the superdeformed band in ^®^Dy. For spher¬ 
ical shape, the two formula for the moment of inertia 
give energy differences (in opposite direction relative to 
the light nuclei) of only around 0.2 MeV for 7 = 60 but 
the difference increases to 0.7 MeV at the calculated 2:1 
deformation of the superdeformed band. 

With the present formalism we are thus able to put ex¬ 
perimental and theoretical results on a common absolute 
scale for high-spin calculations with no pairing. In pre¬ 
vious calculations, either the nuclear binding energy at 
the ground state has been considered as a function of N 
and Z or the energy of a specific nucleus has been consid¬ 
ered as a function of angular momentum, 7. Systematic 
trends or variations of the microscopic energy at high spin 
might provide us with new insight into high-spin phenom¬ 
ena, for example making it possible to compare regions 
where angular momentum is built mainly from collective 
rotation and mainly from single-particle excitations, re¬ 


spectively. A different aspect is that the prediction of 
prompt particle decay from high-spin states, which has 
been observed for example in the nucleus ^®Cu dis¬ 
cussed here, requires reliable estimates of the absolute 
energy difference between the mother and daughter nu¬ 
cleus. 

A natural extension of the present approach is to in¬ 
clude many more nuclei to see if the general trends of Fig. 
IHare still the same. One could also make a new mass fit 
in the high-spin region varying rg and a in the moment 
of inertia formula or in the long run, adjust ‘all parame¬ 
ters’ in a global fit in the N, Z, 7 space in which case it 
would clearly be necessary to also include pairing. One 
possibility might be to calculate the pairing energy only 
for 7 = 0 and postulate a semiempirical formula how it 
slowly disappears with angular momentum in which case 
it would be possible to keep the configuration tracing in 
the high-spin regime which is one of the most important 
features of the CNS approach. Especially for light nu¬ 
clei, the high-spin data appears important. For example, 
it might be possible to obtain a more reliable estimate of 
the stiffness towards deformation for nuclei with A < 100 
where fission barrier data are unknown or uncertain but 
where the energy of several strongly deformed high-spin 
bands are well established. 

In summary, a consistent recipe has been introduced to 
consider the total energy in the full N, Z, I space. Good 
agreement between calculations and experiment is ob¬ 
tained using a standard liquid drop expression with an 
^1/3 curvature term, an average moment of inertia cal¬ 
culated from a diffuse surface mass distribution and shell 
corrections based on the modified oscillator potential. 
For high-spin states in the mass range A = 20 — 200, 
the discrepancies appear comparable to those obtained 
in state of the art mass calculations. The consequences 
of using different models when calculating the rotating 
liquid drop energy were investigated. 
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